Pressure-Driven Orbital Selective Insulator to Metal Transition and Spin State 

Crossover in Cubic CoO 
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The metal-insulator and spin state transitions of CoO under high pressure are studied by using 
density functional theory combined with dynamical mean-field theory. Our calculations predict that 
the metal-insulator transition in CoO is a typical orbital selective Mott transition, where the t2 g 
orbitals of Co 3d shell become metallic firstly around 60 GPa while the e 9 orbitals still remain 
insulating until 170 GPa. Further studies of the spin states of Co 3d shell reveal that the orbital 
selective Mott phase in the intermediate pressure regime is mainly stabilized by the high-spin state 
of the Co 3d shell and the transition from this phase to the full metallic state is driven by the 
high-spin to low-spin transition of the Co 2+ ions. Our results are in good agreement with the most 
recent transport and x-ray emission experiments under high pressure. 

PACS numbers: 71.30.+h, 71.27. +a, 75.30.Wx, 91.60.Gf 
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Although the Mott metal-insulator transition (MIT) 
has been studied extensively for decades, most of the 
works are focused on the single band Hubbard model, 
where the Mott transition is driven completely by the 
ratio of the local Coulomb interaction and band width. 
While most of the Mott MITs in realistic materials [l| in- 
volve more than one band where the transition is driven 
not only by the local Coulomb interaction but also by 
the distribution of the electrons among these bands. For 
example, the redistribution of the four electrons among 
three bands may lead to so-called orbital selective Mott 
transition (OSMT)[1-01 in the t 2g bands. On the other 
hand, in many systems the redistribution of the elec- 
trons among different bands is induced by the crossover 
in spin states, i.e. the high-spin (HS) to low-spin (LS) 
transition [||]. Therefore in realistic materials (e.g. in 
3d transition metal compounds) the Mott MIT and spin 
state crossovers are closely related to each other [9l4l5|. 

Recently the high pressure experiments on charge 
transfer insulator CoO revealed very interesting behav- 
iors in both transport properties and x-ray emission spec- 
troscopy (XES). The transport measurement 16] indi- 
cated that with the increment of pressure there are two 
transitions in resistivity, one happens around 60 GPa and 
the other one takes place around 130 GPa. While the 
room temperature XES measurements 17lll8| on the sim- 
ilar sample show that the spin state of Co 2+ ions persist 
in the HS state all the way to 140 GPa, after which the 
crossover from HS to LS states happens. Hence the in- 
terplay between the HS-LS transition and the two-step 
Mott insulator transition becomes the key factor to un- 
derstand the underlying physics in CoO. 

The pressure-driven MIT and magnetic moment col- 
lapse in transition metal oxides have been studied by 
first principles calculations wi delyfMlEEISl. As 



to CoO, its magnetic state transition under pressure was 
first discussed[13] within the Stoner scenario by employ- 
ing the local spin density approximation (LSDA) and 
generalized gradient approximation (GGA) approaches 
of density functional theory. A transition from HS to 
nonmagnetic metallic state was found around 88 GPa. 
Since these calculations did not take the correlation ef- 
fects of the Co 3d shell into considerations, so the exper- 
imentally observed excitation gap for CoO[18( as large as 
2.5 ~ 2.6 eV was not reproduced completely. Recently, 
Zhang et al. 2l[ reinvestigated the pressure-driven mag- 
netic phase transition in CoO by using LSDA + U ap- 
proach. The HS-LS transition is indeed obtained to be 
°f *2 9 e g ~~ * ^2g e g character, but the electronic structure 
transition is insulator to insulator scheme. This contra- 
dicts with the resistivity data under high pressure fl6jj . 
which shows dramatic change in resistivity indicating the 
insulator to metal transition with pressure. 

In this letter, based on the local density approxima- 
tion (LDA) combined with dynamical mean-field theory 
(DMFT) [22], |23| , we have carried out theoretical calcula- 
tions for cubic phase CoO at different volumes for the first 
time. Our calculations show that CoO under pressure is a 
typical system which has OSMT. At ambient pressure our 
calculation gives correct Mott insulator phase for CoO 
with energy gap being around 2.4 eV. At the first transi- 
tion around 60 GPa the t 2g bands become metallic while 
the e g bands still remains insulating until the pressure 
reaches 170 GPa. Therefore in CoO the exotic orbital 
selective Mott phase (OSMP) with metallic t 2g and insu- 
lating e g bands is stable in a quite large pressure window 
between 60 and 170 GPa, Our LDA+DMFT calculations 
also find that the Co 2+ ions remain in HS state during the 
first transition and the crossover to the LS state starts 
only after the second transition, which is in good agree- 
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ment with both the resistivity 16[ and XES data 171 Il8| 
in CoO. 

The LDA+DMFT calculations 0-12^ in the present 
paper have been carried out by using the pseudopo- 
tential plane-wave method with the quantum espresso 
package 27J for the LDA part and continuous-time quan- 
tum Monte Carlo (CTQMC)[I1 HI as the impurity 
solver for the DMFT part. The single particle LDA 
Hamiltonian is obtained by applying a projection onto 
atomic-centered symmetry-constrained Wannier function 
(WF) orbitals including all the Co 3d and oxygen 2p or- 
bitals, which is described in details in Refliil. That would 
correspond to a 8 x 8 p—d Hamiltonian which is a minimal 
model required for a correct description of the electronic 
structure of CoO due to its charge transfer nature [l|. 

The LDA+DMFT calculations presented below have 
been done for crystal volumes corresponding to values of 
pressure up to 280 GPa. For simplicity, all first prin- 
ciples calculations were performed in nonmagnetic con- 
figuration for rocksalt-type crystal structure with lattice 
constant scaled to give a volume corresponding to applied 
pressure. The Coulomb interaction is taken into consider- 
ations merely among Co 3c? orbitals. In the present work, 
we choose U = 8.0 and J = 0.9 eV, which are close to pre- 
30] . We adopt the scheme proposed 
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in reference [11| to deal with the double counting energy. 
The effective impurity problem for the DMFT was solved 
by the CTQMC quantum impurity solver (hybridization 
expansion version) supplemented with recently developed 
orthogonal polynomial representation algorithm [3lj . The 
maximum entropy method|32j| was used to perform ana- 
lytical continuation to obtain the impurity spectral func- 
tion of Co 3d states. Calculations for all crystal volumes 
were preformed in paramagnetic state at the temperature 
of 290 K. 

In FigUl the evolution of single particle spectrum for 
Co 3d states upon compression is shown. The momentum 
integrated spectral function A{uj) under ambient pressure 
shows well defined insulating behavior for all 3d orbitals. 
However the energy gap for e g orbitals is slightly higher 
than that for t2 g states indicating that the latter orbitals 
are closer to MIT than the former ones. At zero pressure, 
the calculated gap for i 2g orbitals of about 2.1 eV and for 
e g orbitals of about 2.3 eVjWhich agree well with optical 
experimental value 2.6 eV[I||. The LDA+DMFT calcu- 
lations made for small volume values corresponding to 
high pressures gave metallic state for CoO starting from 
60 GPa in agreement with room temperature resistivity 
data[l6j. One can see that ti g orbitals become metallic 
whereas e g ones still remain insulating. This behavior re- 
minds the OSMT scenario as discovered in ruthenatesQ 
at first. The spectral functions for t2 g orbitals in FigQ] 
for pressure values larger than 60 GPa become typical for 
strongly correlated metal close to MIT: well pronounced 
Hubbard bands and narrow quasiparticle peek. However, 
the A(u>) for e g bands remains insulating with Hubbard 
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FIG. 1. (Color online) Single particle spectral function of 
Co 3d states vs pressure obtained in LDA+DMFT calcula- 
tions at room temperature. The spectral function is obtained 
from imaginary-time green's function G(r) by using maxi- 
mum entropy method 32] , and the results are cross-checked 
by usin g re cently developed stochastic analytical continuation 
method 33]. The available x-ray photoelectron spectroscopy 
(XPS) and bremsstrahlung isochromat spectroscopy (BIS) ex- 
perimental data[l|| are drawn in this figure as a comparison. 
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FIG. 2. (Color online) Quasiparticle weight Z of Co 3d states 
as a funtion of pressure. The transition zones are highlighted 
by pink vertical bars. 



bands only but the gap size is strongly reduced compar- 
ing with that of ambient pressure. When the pressure 
exceeds about 170 GPa, the e g states undergo a insula- 
tor to metal transition. As is seen in FigQ] at 187 GPa 
the Mott gaps for both 3d orbitals disappear finally. 

In order to reveal the nature of OSMT in Co 3d or- 
bitals upon compression, we make further estimation 
for their quasiparticle weights by using the well-known 
equation 22j : Z^ 1 = 1 — ^jReE(ia;)| w= o, where ReS(iw) 
is the real part of impurity self-energy function at real fre- 
quency axis. In Fig(5]the calculated quasiparticle weights 



for ti g and e g states as a function of pressure are shown. 
It is apparent that the phase diagram can be splitted into 
three different zones: (1) GPa < P < 60 GPa. The 
quasiparticle weights for both the e g and t2 g orbitals ap- 
proach zero, and the system exhibits insulating behavior. 
(2) 60 GPa < P < 170 GPa. The quasiparticle weights 
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FIG. 3. (Color online) (a) The orbital occupancy of Co 3d 
states as a function of external pressure, (b) The principal 
pressure-dependent atomic state probability of Co 3d states 
obtained by LDA+DMFT calculations. 



for tig states become considerable, whereas those for e g 
states remain very tiny. At this range of pressure, it 
gives an exotic OSMP. (3) P > 170 GPa. The quasipar- 
ticle weights for both t 2g and e g states show a dramatic 
increment with pressure, and the system goes into the 
fully metallic state. 

In Figj3][a) we show the evolution of Co 3d occupan- 
cies and the atomic state probability under compression. 
Due to the charge transfer from O 2p orbitals to Co 3d 
orbitals, the total 3d states occupation number is ~ 7.2, 
which is slightly larger than the nominal value 7.0. At 
ambient pressure, the occupation numbers for t2 g and 
e g orbitals are n(e g ) — 0.54 and nfag) — 0.84, respec- 
tively. Those numbers agree very well with the HS state 
of Co 2+ ion in cubic crystal field with two electrons in e g 
states and five electrons in t 2g states. Over the pressure 
range from to 170 GPa, the occupation numbers for 3d 
orbitals basically remains unchanged. While when the 
pressure is larger than 170 GPa, the occupation numbers 
for ti g and e g states show a dramatic change. The nfog) 
increases to 1.0 and n(e g ) decreases to 0.25 eventually, 
which agree well with the LS configuration of Co 2+ ion 
(t% e l character). Thus the evolution of Co 3d occupan- 
cies with respect to external pressure provides an strong 
evidence for the HS-LS spin state crossover in CoO. 

During the Monte Carlo simulation, we keep track of 
the different atomic state configurations visited and draw 
them as histograms, which give complementary informa- 
tion to the variations of occupancies and magnetic states. 
In Figj3jb), we show the histograms for several upper- 
most N — 7 and N = 8 atomic state configurations. 
It is apparent that at ambient pressure, the HS state 
{t2g e % & = 3/2) makes predominant contribution and 
the contribution from LS state (^2 9 e g> & = V^) and in- 
termediate spin (IS) state (i^e 2 , S = 1) can be ignored 
reasonably. As increasing the pressure, the contributions 
from HS state have dropped and LS state and IS state 
tend to grow, and the total spin magnetic moment will 
decrease as well (see FigfJ|. We note that except for 



FIG. 4. (Color online) Spin state transition in CoO. (a) Evo- 
lution of the local spin magnetic moment of Co 3d states. Here 
M e means effective local moment which can be calculated by 
VTxioc and M t (illustrated by triangle) denotes the previ- 
ous theoretical results 20| . (b) Spin-spin correlation function 
Xioc{t) = (S z (0)Sz{t)} under various external pressures. 



2 configuration, the contribution from another HS 
configuration t\ g e 3 g (S = 3/2) is very little. Thus it is 

confirmed that the HS-LS spin state transition for CoO 
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has a t\ g e g 



1 character, which is consistent with 
In addition, the contribu- 



previous assumptions 17J, [21 1 
tion to magnetic collapse in CoO from HS-IS spin state 



transition {t\ g t 



character) can not be neglected 



as well. 

Next, we concentrate our attentions to the magnetic 
properties of cubic CoO under pressure. In FigfJJa) 
it shows the evolution of the local spin magnetic mo- 
ment with pressure. We note that the effective local mo- 
ment M e is defined through the local spin susceptibility 
y/Txioc, where xioc is defined as xioc = Jq drxiodr) = 
dr(S z (0)S z {T)). As is clearly seen in FigfJJa), under 
compression the local moment decreases slightly from its 
ambient pressure HS value down to about 170 GPa. Fur- 
ther compression rapidly degrades the moment, which is 
accompanied by the redistribution of electrons e g — > ti g 
within the Co 3d shell (see Fig(3]). It should be pointed 
out that at low pressure region our results coincide with 
previously published theoretical results 13 , 13| • 

In FigQJb), the spin-spin correlation function XIoc(t) 
under various external pressure values are illustrated. 
For high pressure, a Fermi liquid phase (LS state) is 
identified. While for low pressure, CoO exhibits a well- 
defined frozen moment phase (HS state). In the Fermi 
liquid phase, the spin-spin correlation function behaves 
as x(r) ~ (T/ sin(TY7r)) 2 for times r sufficiently far from 
t = or /3, respectively. For P > 170 GPa, it displays sig- 
nificant Fermi liquid behaviors, which is consistent with 
the obtained OSMT phase diagram (see Figj2]) for CoO. 
The frozen moment phase is characterized by a spin-spin 
correlation function that approaches a non-zero constants 
at large times, as is seen from to 170 GPa. Thus this fig- 
ure reveals a phase transition between a LS Fermi liquid 
metallic phase and a HS Mott phase with frozen moments 
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again. 

Now we discuss the most interesting issue of this paper, 
the relationship between OSMT and HS-LS transition. 
The detail analysis of the spin state indicates that the 
HS-LS transition is the most important driving force for 
the OSMT in CoO. Once the Co 2+ ions keep in HS state, 
the e g bands are always half filled, which greatly favors 
Mott insulator phase and very hard to become metallic. 
On the other hand, in the HS state the filling factor of 
t2 g bands is only 1 /6 (in terms of hole density) , which 
makes it much easier to become metallic. For example, 
the previous studies on multiband Hubbard model [34] 
show that the critical U c for half-filled two-bands model 
is around 1.7W and that of the 1/6 filling three-bands 
model is around 2.1W, where W is the band width. Thus 
apparently it is the different situation in filling factors 
between e g and t2 g bands which makes the two sub-shells 
behave so differently under pressure. Once the LS state 
stabilized above 170 GPa, the e g orbitals are no longer 
half filling and turn to the metallic phase eventually. 

In summary, we conclude that the two-step like Mott 
transition in CoO can be understood as a typical OSMT. 
At ambient pressure CoO is a typical Mott insulator with 
energy gap around 2.4 eV. At the first transition around 
60 GPa the t?. g bands become metallic while the e g bands 
still remains insulating until the pressure reaches 170 
GPa. Therefore in CoO the intriguing OSMP with metal- 
lic ti g and insulating e g bands is stable in a quite large 
pressure window between 60 and 170 GPa. Our theo- 
retical calculations also find that the Co 2+ ions remain 
in HS states during the first transition and the crossover 
to the LS states starts only after the second transition, 
which is in good a gree ment with both the resistivity [16| 
and XES data[I3,I3 for Co °- Further analysis of the 
calculated results show that the HS-LS transition is the 
main driving force of the Mott MIT in cubic CoO. 
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